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fvq Abstract 

-.^ A theoretical study has been carried out to analyze the available results from the inelastic 

fvj neutron scattering experiment performed on a quasi-two dimensional spin-^ ferromagnetic material 

K2C'uF4. Our formalism is based on a conventional semi-classical like treatment involving a model 
of an ideal gas of vortices/anti- vortices corresponding to an anisotropic XY Heisenberg ferromagnet 



Oh 



Xy^ on a square lattice. The results for dynamical structure functions for our model corresponding to 

^_ spin- 1 , show occurrence of negative values in a large range of energy transfer even encompassing 

I the experimental range, when convoluted with a realistic spectral window function. This result 

indicates failure of the conventional theoretical framework to be applicable to the experimental 

'~~ ' situation corresponding to low spin systems. A full quantum formalism seems essential for treating 

D such systems. 

;_( 

j/3 PACS: 78.70.Nx-inelastic neutron scattering in condensed matter, 75.10.Jni"Heisenberg model, 

-<— » 75.30.Kz-Kosterlitz-Thouless transition in magnetic systems 

S ... 

_L 1 Introduction 

O Low dimensional and in particular two dimensional magnetism has attracted a great deal of interest 

I ^1 in the past three decades [HEllS]. In particular, in one dimension the existence of both solitonic and 
spin wave excitations were thoroughly studied through inelastic neutron scattering experiments as 
\^ well as theoretical analysis for CsNiF^ ^. Similar studies were carried out searching for topological 
Q\ excitations in various quasi-one dimensional systems which are almost ideal realization of nearest 
^O neighbour Heisenberg antiferromagnetic chain [5j. 

cn 

^J In many of the above systems the experiments showed central peak (peak corresponding to a; = 0) 

(^ in the dynamical structure function when plotted in constant "g" scan. This motivated the ex- 
CN perimentalists further to investigate two dimensional and quasi-two dimensional magnetic materi- 

. . als. With the availability of improved quasi two-dimensional ferromagnetic and anti-ferromagnetic 
.^ materials investigations along this line has become todays one of the primary interests both theo- 
S^ retical and experimental. These include layered systems such as K2CUF4, Rb2CrCl/i, magnetically 
^ intercalated graphites such as C0CI2, layered ruthenates, layered manganites and high Tc cuprates 
[BISlEllllinilTllHllHl dni [m [12]. Moreover large amount of information on the spin dynamics, ex- 
tracted from inelastic neutron scattering are available. Advances in numerical and computational 
techniques have also contributed to the understanding of both spin wave and topological excitations 

[laiiiiiiaiisiiiT]. 
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On two dimensional magnetic systems the concept of topological order was introduced by Koster- 
litz and Thouless and independently by Berezinskii [18^ [19] . Their ideas back-ed by analytical and 
numerical calculations led to the proposal for the existence of topological vortices and anti-vortices in 
a typical ferromagnetic XY model on a two dimensional lattice. According to these ideas vortices and 
anti-vortices are frozen as bound pairs below certain transition temperature called Tkt or Tbkt and 
above this temperature they become mobile and nearly free |18[ [19] . 

In this work we initiate a theoretical investigation regarding the applicability of a semi-classical like 
treatment of the dynamics of topological excitations to the inelastic neutron scattering results for real 
systems [6l [T3[ [T^ [20j . In recent years inelastic neutron scattering experiments are being done mostly 
on layered ruthenates like Ca2-xSrxRu04, layered anti-ferromagnet like CuGeOs, layered manganites 
and some layered cuprates [D [2 |3l [lOl [TTl [21]. The layered anti-ferromagnet, CuGeO^ being spin- 
Peierls compound, is proposed to exhibit a Berezinskii-Kosterlitz-Thouless transition in the vicinity of 
spin-Peierls transition temperature. The two dimensional spin half XY model was investigated and the 
validity of Berezinskii-Kosterlitz-Thouless transition was confirmed |22j . The existence of Berezinskii- 
Kosterlitz- Thouless transition was proposed long ago in K2CuF4{S= ^ layered ferromagnet) [6]. 
According to our knowledge layered ferromagnets with spin- 2 are the least studied systems, both 
from theoretical and experimental point of view, till date. 

An extensive experimental study of spin-dynamics in a layered ferromagnet has been carried out 
by Hirakawa et. al. [6] using neutron scattering probe on K2CuF4^. Their results exhibit a central 
peak (at a; = 0) in the plot of "neutron count vs. frequency" at a fixed value of the wave-vector 
q. Subsequent developments of approximate analytical theories and Monte Carlo Molecular Dynam- 
ics(MCMD) analysis have suggested that the existence of central peaks is partly due to scattering of 
neutrons from moving vortices and anti- vortices |14[ [T5] . 

Here we aim to examine how far the picture of ideal gas of vortices and anti-vortices could be 
extended to the quantum spin models. For this purpose, we choose K2CUF4 as the reference system. 
This is a spin- 2 quasi-two-dimensional ferromagnetic material, on which extensive neutron scattering 
studies have been done. 

The plan of the paper is as follows. In Sec. 2, we briefly describe the classical theory of mobile 
vortices and anti-vortices. In the same Section we explain our mathematical formulations in detail. 
In Sec. 3 we discuss our calculations and results. In Sec. 4 we present the conclusions and the future 
plan. 

2 Mathematical Formulation 

The dynamics of mobile vortices in a ferromagnetic system has already been treated both analytically 
and numerically by Huber [1^ and Mertens et. al. [2]. In this work we apply their classical 
formalism to study the phase transition in K2CUF4. We calculate the spin-spin correlations taking the 
experimental situations into account. For our purpose we present a brief description of the analytical 
treatment developed in Ref. [14j. The starting Hamiltonian is, 

n = -jJ2{sfs^ + sfs] + xsfs^), (1) 

where i,j label the nearest neighbour sites on a two dimensional square lattice, J is the coupling con- 
stant and the classical spin vector is Sj = {Sf, Sf , Sf). This is an anisotropic Heisenberg Hamiltonian 
which, for J > 0, represents a ferromagnetic system. The quantity A is the anisotropy parameter 



whose XY and isotropic Heisenberg limit correspond to A = and 1 respectively. The general time 
dependent spin configuration in spherical polar coordinate system is given by, 

Sx = S cos(l){r,t)sin6{r,t), 
Sy = S sincj}{r,t)sin6{r,t), 
S, = Scose{r,t), (2) 

with r = {x,y). Following the formulation of Hikami and Tsuneto, the solutions are given by |20j . 

(j) = ibarctan(|) and 

e = '"(lie"'^/^") for r^ry, 

= or vr r — ;• 0, (3) 

for single vortex centred at r = (0, 0), where (3) describes the asymptotic behaviour of 6. Here vortex 
core radius is gi^ 

the spin vortex. 

The definition of the spin-spin correlation function is given by, 



core radius is given by [15], r^ = , " This type of spin configuration defines a 'meronic' type of 

-y/2(l-A) 



5(r,t) = (S(r,t) -8(0,0)) 

= (5-(r,t)5-(0,0)) + {Sy{v,t)Sy{Q,Q)) + (5^(r,t)5^(0,0)), (4) 

where (...) represents the thermal average. In the case of classical gas of ideal vortices the ther- 
mal average has to be done by taking Maxwellian velocity distribution function. Here S^^{r^t) = 
{Sx{r,t)Sx{0,0)) and 5^^(r, t) = {Sy{r,t)Sy (0,0)) are in-plane correlations and 5^^(r, t) = {Sz{r,t) 
Sz{0,0)) is the out-of-plane correlation. 

The effective analytical expression for the in-plane correlation can be taken as |14j . 



S--(r,t) = ^exp{ 



2 



2.2 



2+rt 



ni/2 

}, (5) 



with 7 = ^^^, where u is the root mean square velocity. Here ^ = Co^'' is the vortex- vortex 
correlation length. The root-mean squared velocity of the vortices was first calculated by Huber as 

where n^ is the density of free vortices at T > Tkt- The Fourier transform of S^^{r,t) in (5) gives 
rise to the in-plane dynamical structure function given by, 

^"""^ (q, ^) = ^ [a;2 + ^2(l + ^2g2)]2 • (7) 

This is a squared Lorentzian, peaked at w = 0, with q dependent width, 

r = ^{vr(V2-l)}V2(^iyrTev). (8) 

Exactly same results holds for S'^^(q, w) also. 



From the definition of S'^^(r, t), it can be sliown tliat the out-of-plane correlation is given by [14j . 
S^'(r,t) = n{,S^ f f d^Rd^uP{u)cos9{r-K-ut)cos9{R), (9) 



where P{u) is the Maxwell velocity distribution for a single vortex. Performing first the spatial Fourier 
transform and then the temporal Fourier transform, it can be shown that, the out-of-plane dynamical 
structure function has the form, 

Here |/(q)| is the velocity independent vortex form factor and it has the form /(q) = J (Pr cos 6{r)e~^^''' . 
The form of ^^^(q,!^) as in eg"'(10) exhibits a central peak at u; = 0. The width of the central peak 
is r^ = uq i.e., linear in q. 

In a typical inelastic neutron scattering experiment the count rate is related to the dynamical 
structure function as [23], 

/(q, w) oc ^(q,^). (11) 

At finite temperature there always exist creation and annihilation of excitations. A detailed balance 
condition is always needed to relate the intensities of up scattering (huj < 0) and down scattering 
{huj > 0). True quantum mechanical 5'(q, w), denoted by 5z)B(q, w) is recovered by the relation , 

Where the factor ^ .^^ is called the Windsor factor [T7]. This Sdb{^, ^) incorporates the detailed 

balance condition, as required by the thermal equillibrium. Another important factor, which has to 
be taken into account, is the instrumental resolution function R{t) or R{uj — u'). This essentially 
incorporates the different independent instrumental properties that affect the incident and scattered 
beam of neutrons f23] . In order to compare theory with experiment one has to convolute the theoretical 
expression, obtained from a model under consideration, by the resolution function |24j . Thus we 
consider the convoluted dynamical structure function Sconv. (q, w) given by, 

SconvXci, ij) = J dt J d\R{t)S{r,t)e'^^-''-'^'^ 

= R{q,uj — uj')S{q^,uj')doj' . (13) 

The resolution function has to be chosen so as to give minimum ripples at the end points of the 
resolution width. For this purpose a suitable window function such as Tukey window function may be 
chosen, 

R{t) = - [1 + cos(2^t/t^)] for \t\ < tm/2 

= otherwise (14) 

The parameter tm, occurring in the Window function, can be set from the resolution half width 
obtained from experimental data. 



2.1 In-Plane Dynamical Structure Function 

In our formulation for the in-plane dynamical structure function we take into account the Tukey 
window function, as mentioned above (see equations (14)). Using (5), (13) and (14) we compute the 
Fourier transform of in-plane spin-spin correlation, 

ei(q-r-^t)^ (15) 

Now, J cPre'^'^'^ = J^ r dr J^^ dOe^'^^'^"'^^ = /(^ r dr Jo{qr), where Jo{qr) is Bessel function of order 
zero. The spatial integration is performed from zero to a certain radius Rq. A final expression for the 
convoluted in-plane dynamical structure function takes the form, 

rJo{qr)R{t)cos{ojt). (16) 

Since, S^^{Y,t) and R{t) are both even function in t , only cos{ujt) contributes to the temporal part of 
the integration. From symmetry Y component of the in-plane dynamical structure function S}!^^^{(1, to) 
is same as X component of the in-plane dynamical structure function S'^^^(q, tj). Let us note that 
in the above analysis the formulation holds only for T > Tkt- For T < Tkt the vortex- vortex 
correlation length ^ is not defined and hence the formalism can't be extrapolated below Tkt- 

2.2 Out-of-Plane Dynamical Structure Function 

The out-of-plane dynamical structure function is given by, 

S'conM, OJ) = j R{U- J)S^\c^,J)dJ , (17) 

where i?(a; — uj') is the Fourier transform of R(t). The reason for taking (17) as the expression for 
convoluted out-of-plane dynamical structure function is that, unlike (5), an analytical expression for 
^^^(r, t) can't be evaluated from (9). So one has to start from eg"(10). 

The integral in (17) has been computed numerically. The S'^o„^(q, oj) defined here corresponds only 
to the mobile vortices; whereas the experimental data contain the contributions from bound vortices 
and fragile 'spin wave like' modes. These fragile 'spin wave like' modes are the largely decaying spin 
wave modes above the ferromagnetic- paramagnetic transition temperature (Curie temperature). In 
order to compare with the experimental observations, one has to extract the mobile vortex contribution 
from the experimental data. This can be done by subtracting the fragile mode contribution and 
the frozen vortex contribution from the experimental data. The fragile mode contribution has been 
subtracted by taking the fragile mode contribution above transition temperature to be same as the 
spin wave contribution just below transition temperature. This is valid as long as we are considering 
the temperature which are not far below or above from the transition temperature. To find the 
approximate analytical expression for S^^{(\,uj) due to bound vortex contribution, the limiting value 
of u is taken as n — t- in (10). Then from (10) it easy to find an expression for 5^|^„„^(q, oj) namely, 

S-f,^U^,u:) = ^nl\f{q)\H{u), (18) 

Where n^ is the bound vortex density. Since, the system has no net topological charge we can assume 
that there are equal number of vortices and anti-vortices present in the system and we can take 



nl + Tiy = ^ assuming square lattice structure. This is correct as long as the temperature is just 
below Tkt where all the vortices are frozen but once the temperature crosses Tkt some of the bound 
vortices become mobile and the bound vortex density can be approximated as, 

nl^i\-nf), (19) 

where, n^ is in the units of inverse of plaquette size(a^). Since, nl ~ ^q exp(— 26/-y/T) [H], n^ given 
by (19) is temperature dependent. Here ^o is of the order of lattice parameter. Using (18) and (19) 
the bound vortex contribution has to be subtracted carefully from the experimentally observed count. 

We would like to point out that we could not apply this procedure to extract out bound vortex 
contributions in the case of in-plane dynamical structure function (see Sec: 3). 

2.3 Total Dynamical Structure Function (Spin-Spin Correlation) 

The general expression for the total dynamical structure function is, 

5(q,^) = ^-L.y'dVy'cit 5(r,i)i?(t)e^('i--*). (20) 

where, the total spin-spin correlation isS'(r, t) is defined by (4). So, the total dynamical structure 
function is 5(q, w) = 5^^(q, cj) -|- ^^^(q, a;) -|- ^^^(q, w). Since, X and Y components of the spins are 
symmetric we have, S^^{c[,uj) = ^^^(q, Ci;)and the total dynamical structure function takes the form, 

5(q,^) = 25""(q,a;) + 5^^(q,w). (21) 

Here, we would consider (21) only for mobile vortices. 

It is an important fact that the formalism explained above incorporates the Windsor factor and the 
presence of h in the quantum expression of magnetic moment corresponding to the spins constituting 
the vortex [13]. Therefore the formalism looks like a semi-classical one. Henceforth we will call our 
combined theoretical approach 'semi-classical like'. 

3 Calculations and Results 

We apply the formalism of Sec: 2 on a real material K2CUF/1 for which neutron scattering experiments 
have been performed [^. It is a quasi-two-dimensional spin- 2 ferromagnet, where the interaction is 
mainly Heisenberg type with only 1% X-Y like anisotropy . The transition is close to KT type with 
slight modification due to Heisenberg type interaction. The magnetic lattice structure for K2CUF4 
is approximately a body centred tetragonal lattice i.e. a lattice, composed of stacking of 2D square 
lattices [6]. The physical parameters are given in the Table 1, which have been used throughout the 
calculation. 

We start with the investigation of the in-plane correlation (in-plane dynamical structure function) 
5^^(q,a;). The radius Rq in (16) is {VmFTmfi)a for a 100 x 100 lattice, as used in the MCMD 
analysis by Mertens et. al. [H], where a is the lattice parameter. We set the value of tm according 
to the experimental resolution width (0.01 meV) [6]. We compute ^^^^^(q, cj) numerically, for two 
different temperatures, 6.25 K and 6.75 K, for q(planar)= 0.04 reciprocal lattice units(in the units of 
^), experimentally fiq being the momentum transfer. There are two threshold values of g [6], namely 
qi = 0.06 and g2 = 0.01, where for q > qi the system behaves like 2D Heisenberg system and for q < q2 
the system behaves as 3D XY system. For q2 < Q < Qi the system behaves as 2D XY system. We have 



Table 1: Relevant parameters for K2CUF4 [H] 



parameter 



magnitude 



exchange coupling (J) 11.93 K 

lattice parameter (a) 4.123 A 

'b' 1.5 

'Tkt 5.5 K 



T= 6.25 K, q= 0.04 




-0.2 

energy transfer ho in meV 



0.4 



Figure 1: Comparison between the convoluted in-plane dynamical structure function 5^^j,(q, w) (eqn. 
16) and unconvoluted in-plane dynamical structure function 5^^(q, w) (eqn. 7)at T = 6.25 K and q = 
0.04. Solid line is for convoluted theoretical expression and dotted line is for unconvoluted theoretical 
expression(squared Lorentzian). ^ = 58.09a, u = 0.0614^-^—, and width T^x = 0.0012?7T-ey for squared 
Lorentzian at T= 6.25 K. 



T= 6.75 K, q= 0.04 




-0.2 

energy transter ho in meV 



0.4 



Figure 2: Comparison between the convoluted in-plane dynamical structure function 5^o„j;(q, w) (eqn. 
16) and unconvoluted in-plane dynamical structure function S^^{c[,uj) (eqn. 7)at T = 6.75 K and q = 
0.04. Solid line is for convoluted theoretical expression and dotted line is for unconvoluted theoretical 



expression(squared Lorentzian). ^ 
Lorentzian at T= 6.25 K. 



22.25a, u = 0.1352^ and width F^^ = 0.0035mel/ for squared 



varied the energy transfer hu from -0.3 meV to +0.3 meV, which includes the range -0.2 meV to +0.2 
meV as taken in experiment jBj. The convoluted in-plane dynamical structure function is plotted in 
Fig. 1 and Fig. 2, where tnat = js(s+i) ^® ^^^ natural time unit for the system/material (in our case 
K2CuFii). These figures indicate that after convoluting with the Tukey window function, the in-plane 
dynamical structure function no longer remains squared Lorentzian, though in both the cases central 
peaks persist. The width of the S^^^^{c{,oj) curve is much larger than that of the squared Lorentzian. 



We notice that the convoluted in-plane dynamical structure function function 5|f^„^(q, oj) has become 
negative just above 0.1 meV. The occurence of negative values of the dynamical structure function 
has been dealt with in detail in sec. -4 and in Appendix. 



Again comparing Fig. 1 and Fig. 2 we find that the width of the squared Lorentzian increases with 
the increase of temperature whereas that of the S'^^^(q, w) does not undergo any change. Later, we 
will present a comparison of the convoluted total dynamical structure function with the experimental 
one(See Fig. 5 and Fig. 6). 



We now evaluate the out-of-plane dynamical structure function ^^^^^(q,^) for two different tem- 
peratures, 6.25 K and 6.75 K, for q(planar ) = 0.04 r.l.u, using (17). The expression for Riu — u') 

is, 



R{oj 



ijj 



1 

47r' 



-sm 



{UJ - J)tr 



1 



1 



UJ 



UJ UJ 



(22) 



Uj' + 27: /tm UJ — Uj' — 2tX jtm 

(q, o;) . Here also the reasons for the choice of temperatures 
and q(planar) are same as that for the in-plane dynamical structure function. In Fig. 3 and Fig. 4 



We use the same value of tm as used for S^^^^ 



J^ 1 1 .93 K. q{plaiiar)^ 0.04, T^ 6.25 K 
circles are observed dala wilh bound voiles contribution subtracted, solid line is theoretical plot 



400 



300 - 



200 



■^ >, 100 ■ 



-100 



-200, 



— t> — 



-0.2 



-0.1 



0.1 



0.2 



h(0 in mcV 



Figure 3: circles are observed (experimental) data, where contributions from the fragile modes as well 
as the bound vortex contributions have been subtracted & solid line is the plot of properly convoluted 
out-of-plane dynamical structure function ^^^(qja;) (theoretical). ^ = 58.09a, u = 0.0614^r^. 



we have plotted the out-of-plane correlation, 5^o„„(q, w). We have varied the uj' from 
(17), where tm is estimated from the resolution width as before. 



to 



m 



In this case the bound vortex contribution has been subtracted carefully, using (18) and (19), from 
the observed count at 6.25 K to obtain the effective mobile vortex contribution. The methodology for 
extracting the mobile vortex contributions from the experimental data has been explained in Sec 2.2. 
As long as the counts at 6.75 K are concerned, the fragile 'spin wave like' modes are highly decaying 
so that it can't be assumed to be the same as the true spin wave modes observed at 5 K. So only 
bound vortex contribution has been subtracted at 6.75 K. The normalization factors, required for the 
quantitative comparison between the theoretical and the experimental results, have been estimated 
from the neutron count extracted from the experiment on K2CuF4^ j6]. 

We find that the out-of-plane dynamical structure function is also negative within the resolution 
width(see Figs. 3 and 4)! Moreover experimental peak is out side the resolution width, while the peak 
corresponding to the S^oni;(^i^) is at cj = 0. 

The above calculations lead us to the theoretical estimate for the convoluted total dynamical struc- 
ture function 5*°^'!j|(q, w) given by (21). In Figs. 5 and 6, S'*°^'^'(q, w) has been compared with the 
filtered experimental data obtained by subtracting the bound vortex contributions and fragile 'spin 
wave like' contributions (see Sec. 2.1 and 2.2). In these plots the intensities of the experimental peak 
and that of the central peak of the 5*°^'*^'(q, w) have been matched. 



It is clear from Fig. 5 that at 6.25 K the experimental peak occurs approximately at 0.08 meV, 
which is way outside the resolution width. At 6.75 K [see Fig. 6] the peak of the experimental graph 
is not far from the central peak. It is reasonable to say that as the temperature is increased, we 



J^ 1 1 .93 K. q{plaiiar)^ 0.04, T^ 6.75 K 
circles ait observed dala with bound vones contiibution subtracted, solid Hue Is theoretical plot 
300 r 



200 



a i 

2 I 'H 



-100 



I n " " ' 



-0.2 



-0.1 



0.1 

Ti(0 in mcV 



0.2 



Figure 4: circles are observed (experimental) data, where only contributions from bound vortices 
have been subtracted & solid line is the plot of properly convoluted out-of-plane dynamical structure 
function 5^^(q,a;) (theoretical). ^ = 22.25a, u = 0.1352^. 
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S 400- 



^ 300 



200- 
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K 100- 



theoretical plots 

• experimental data 
q(pianar) ^ 0.04. T= 6,25 K 




-0.3 -0.2 -0.1 0.0 0.1 0.2 

energy transfer hra/ 2n in meV 



0.3 



Figure 5: total dynamical structure function S'*"^" (q, w) at T = 6.25 K and q = 0.04- solid line is for 
convoluted theoretical results and dots are filtered experimental data. ^ = 58.09a, u = 0.0614^-^. 
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• experimental data 
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Figure 6: total dynamical structure function 5*°^" (q, w) at T = 6.75 K and q = 0.04- solid line is for 
convoluted theoretical results and dots are filtered experimental data. ^ = 22.25a, u = 0.1352t^^. 



are getting better agreement of the S'*°„° (q, w) with the experimental observations. This agreement 
isregarding the position of the central peak. Apart from the central peak there are two other peaks at 
finite frequency at both the temperatures. These are nothing but the reminiscent of the out-of-plane 
dynamical structure function contribution as seen from Figs. 3, 4, 5 and 6. This signifies the fact that 
the in-plane correlation is largely dominating over the out-of-plane correlation. 



The total spin-spin correlation is still negative just above 0.1 meV. Though it is true that the 
dynamical structure function can't be negative, here in our case the negativity occurs as a result of 
the convolution of analytical expression of ^(q, w). Even for a conventional long range ordered system, 
the dynamical structure function corresponding to a classical pure spin wave comes out to be negative 
beyond a certain range of frequency when convoluted with any spectral window function. Furthermore 
the above peculerity persists even when quantum effects are incorporated through a detailed balance 
factor [see Appendix]. 

The inclusion of quantum mechanical detailed balance factor in the semi-classical like treatment 
for dynamics of mobile vortices and anti-vortices, is not even causing any appreciable asymmetry, as 
seen in the theoretical plots in our case of spin-^- The theoretical plots are largely symmetric around 
a; = 0. A very small asymmetry in the theoretical plots are being seen for higher values of u while 
the experimental data are showing clearly the asymmetry. 

It may be noted that in our analysis the bound vortex contributions have been approximately 
estimated only for out-of-plane dynamical structure function ^^^(q, w). This is because, in this case, 
we are able to truncate the expression, as given in (10), to the regime T < Tkt, by making u — )• 0. 
In (10), there exists no explicit dependence of 5^^(q, w) on the correlation length ^. In case of in- 
plane correlation, as given in (5), we need to find ^ for T < Tkt due to its explicit appearance in 
that expression. Since ^ is not defined for T < Tkt we are not able to estimate the bound vortex 
contribution for in-plane dynamical structure function. 
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In summary, we find that the width of the convoluted in-plane dynamical structure function is much 
larger than that of the squared Lorentzian. Values of the in-plane dynamical structure function comes 
out to be negative beyond a finite range of energy transfer. The convoluted out-of-plane dynamical 
structure function becomes negative as well; however this happens within the resolution width about 
the central peak (peak at cj = 0). The total convoluted dynamical structure function also becomes 
negative in the regime where the in-plane dynamical structure function had become negative. No 
appreciable asymmetry is created even after including the Windsor factor. We find that for both the 
temperatures the convoluted total dynamical structure function is symmetric around oj = 0; whereas 
the experimental observation is not. The theoretical model of semi-classical treatment of ideal gas 
of unbound vortices tends to agree with the experimental observations better at higher temperatures 
(for spin-2 system), when we consider the experimental results at T= 6.25 K and T= 6.75 K (Figs. 5 
and 6). It is worthwhile to point out that same results hold for unbound anti- vortices also. 

4 Conclusions &; Discussions 

The laws of quantum mechanics which govern all real systems, ensure the dynamical structure functions 
to be always positive definite [23]. We find in our analysis however, that the semi-classical treatment 
based on ideal gas of vortices (anti- vortices) for a low spin system leads to the occurence of negative 
values of dynamical structure function, over a large range of energy transfer, when convoluted with 
any standerd resolution function. 

Based on the analysis carried out in the Appendix we can infer that for the dynamics of mobile 
vortices and anti- vortices, the negative values of SconvX^^ ^) ^^^ occurring due to the following factors; 
i) the choice of the resolution function, whicn in our case is the Tukey function, 

ii) the choice of the value of resolution width Aoj, which in this case is made fixed by experimentally 
imposed resolution width, 

iii) use of a semi-classical like treatment to extend the classical theory of dynamics of mobile vortices 
(anti-vortices) to a quasi two dimensional spin ^ ferromagnet which is quantum mechanical. 

To avoid the negativity in the Sconv. (<!) ^) we could have chosen a different resolution function. 
Indeed, it has been shown that most of the resolution functions are more or less oscillatory in the 
fourier space |24^ [29] . An extra smoothening factor can be used to dampen the oscillation of the 
resolution function. This extra factor is eventually related to the resolution width and it makes 
the resolution function smoother if the resolution width is decreased [23] . However, in our case the 
resolution width is fixed from the experiment and consequently the oscillation of the resolution function 
can't be avoided by merely changing the resolution function. 

It has been found that there is a range of uj over which Sconv. (Q) ^) remains positive. We can call 
it as the physically admissible range. This range is related to the magnitude of the spin occuring in 
the theoretical model under consideration and to the resolution width. On the basis of the analysis 
(presented in the Appendix) it is expected that even in the case of dynamics of mobile vortices and 
anti-vortices the physically admissible range would be larger for higher spin value and smaller for lower 
spin values. 

Another way to avoid the negativity is to assume, SconvX^-, w) = outside the physically admissible 
range [M]- If this physically admissible range is within the range of experimental interest then the 
assumption is not applicable. In our case of spin-half ferromagnet the physically admissible range is 
well within the range of uo over which the neutron scattering data has been taken in the experiment 
(as seen from Figs. 5 and 6). 
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Hence, the negative values of Sconv.i^?'^) can only be due to the use of the semi-classical like 
treatment to extend the classical theory of dynamics of mobile vortices and anti-vortices to a quasi 
two dimensional spin ^ ferromagnet {K2CUF4). 

Moreover, the convoluted out-of-plane dynamical structure function computed from our semi- 
classical like treatment, becomes negative within the experimentally imposed resolution width itself. 
Thus the central peak occurring in this case may not posses a well defined width. 

The agreement between the behaviour of dynamical structure functions obtained from our theoretical 
calculations and that from the experiment, in terms of the peak position and the overall shape, is found 
to be fairly good at temperatures much larger than Tkt- 

Although the vortices (and anti-vortices) are extended objects the Maxwell- Boltzmann distribution 
can still be used for the motion of the centre of mass of these objects. 

Our investigation presented in this paper brings out the fact that a complete quantum mechanical 
treatment is essential for describing the detailed features of the dynamics of unbound spin vortices and 
anti-vortices corresponding to low spin magnetism systems. As a first step towards this, a theoretical 
framework for describing static quantum spin vortices and anti- vortices and their topological properties 
has been developed [S^ [Ml EZJ [2H] . An extension of this formalism to the case of mobile spin vortices 
and anti-vortices is crucial for the quantum mechanical calculation of dynamical structure function. 
This would go a long way towards an explanation for the experimental results observed for the genuine 
quantum spin systems like K2CUF4. 
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Appendix: Sconv{<^, ^) corresponding to classical spin- wave 

For classical spin wave at very low temperature corresponding to a classical Heisenberg ferromagnet, 
the dynamical structure function has the form:- 

5(q, w) = (5(cj^ - UJg) 

= ^6{io-u;g) + d{uj + ujg)]. (23) 

ZUJq 

Where, ojq = h^^J^S{S + 1) z{l — 7^) for a cubic lattice. Here z is the number of nearest neighbours 
and 7g = ^ErCOs(q-r). 

Convoluted dynamical structure function, Sconvi^i, w) has been defined in (13), where R{u} — a;') is 
the fourier transform of a siutably chosen spectral function. So far Tukey function has been used and 
the same is used in this case also. The fourier transform of Tukey function is given by (22). Using 
(13), (22) and (23) we find the convoluted dynamical structure function corresponding to spin wave, 
S^Zi^, ^) to be, 

5f,^„(q, a;) = ^[fi(t^ - u,) + R{u; + u:,)]. (24) 
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Analyzing (24) we find that there exists a range 

|w| < K + ^) (25) 

within which ^^^(q, lo) remains positive and outside, it becomes negative. This remains true even if 
we choose any other resolution function. This is basically due to the fact that when fourier transform 
is performed on spectral functions (defined in time domain), the resulting functions in (q, w) space 
mostly turn out to be oscillatory |29| . 



One way to avoid these negative values of Sconvi^i w) would be to assume Sconv{^, w) = outside 
If^l ^ i^q + p") |211- This prescription however, can't be taken into consideration when comparing 
theoretical predictions with experimental results, if the energy range of experimental interest contains 
the above mentioned range of w. Another way to avoid the negative values of Sconvi^i, ^) is to decrease 
the value of t^ [see (25)], where tm is related to the experimental resolution width (Acj) by the relation, 

Acj being in energy units. Decrease in tm means increase in the value of Alo which signifies a poor 
experimental resolution width. In comparing the theoretical predictions with the experimental results, 
tm is determined by the relation (26), where Aw is fixed from the experiment. 

Another interesting feature of (25) is the fact that the region of a; where Sconvi^, w) remains positive, 
depends on the spin value S via the relation, 

H < {h'^J^JsiS + 1) z[l - 7,] + ^). (27) 

tm 

This shows that when tm is fixed from the experimental resolution width, S remains as the free 
parameter to determine the range of u as mentioned in (25). Increase in S will increase the range of u 
where Sconv{<i, oj) remains positive. Hence for a system with high spin value S, the range of positive 
Sconv{fl, w) may actually be at par with the energy range of experimental interest, unlike the case of 
spin 2. 

It can be shown from (24) that the range of w, given by (25), will not change even if a detailed 
balance factor (in this case Windsor factor) is introduced to incorporate the quantum effects. 
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